AD-763  717 


r 


IMAGING  OF  EXTENDED  OBJECTS  THROUGH 
A  TURBULENT  ATMOSPHERE 

Hal  T.  Yura,  et  al 


Aerospace  Corporation 


Prepared  for: 

Space  and  Missile  Systems  Organization 
28  June  1973 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


AD  7  63717 

'  •  -  ’’  rs 


a  ,  -  -  AEROSPACE  REPORT  NO- 

AIR  FORCE,  REPORT- NO.  TR-0073(3230-30)-4 

SANSO-TR-73:213  -  - 


Imaging  of  Extended  Objects  Through  a 
Turbulent  Atmosphere 


'Preparcdby-R.  F.  lUTOMlRSKl 
Pacific  Sierra  Research  Corporation 
ami 

H.  T-  YURA' 

Electronics  Research  Laboratory 
Thc'Acrospacc  Corporation 


Prepared  for -SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
AIR  FORCE  SYSTEMS. COMMAND 
LOS  aMGELES  AIR  FORCE  STATION. 

Los  Angeles,  California 

Stprodutf O  bf 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  S  E«DO'*',*<r  1  of  Co^^erce 
S?  t‘<i  v  A  2215! 


APPROVED  FOR  PUBLIC  RELEASE.  DISTRIBUTION  UNLIMITED 


LABORATORY  OPERATIONS 


The  Laboratory  Operation*  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigation*  necessary  for  the  evaluation  and 
application  of  scientific  advance*  to  new  military  concept*  and  system*.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  m  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aeroohysic*  Laboratory.  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-pouer  gas  Lasers. 

Chemistry  and  Physics  Laboratory  Atmospheric  reactions  and  atmos- 
phenc  optics,  eKerrucalreactioa*  w  pofiLted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surfa*  e  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics,  communication  sciences,  applied  electronic*,  semi¬ 
conducting.  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging,  atmospheric  pollution,  millimeter  wave  and  far-infrared  technology 

Materials  Sciences  Laboratory  Development  of  new  materials,  metal 
matrix  composite*  and  new  form#  of  carbon,  test  and  evaluation  of  graphite 
and  ceramic*  in  reentry,  spacecraft  materials  and  electronic  component*  in 
nuclear  weapons  environment,  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue-induced  fractures  m  structural  metals. 

Space  Physics  Laboratory  Atmospheric  ard  ionospheric  physics,  radia- 
t.on  from  the” itmospherc,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow.  magnetosphenc  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere  solar  physics,  studies  of  solar  magnetic 
field*,  space  astronomy,  x-ray  astronom>,  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  spac^  systems 

THE  AEROSPACE  CORPORATION 
El  Segundo,  California 


MSCLilHEl  NOTICE 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


UNCLASSIFIED 


i _ Security  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 

(Security  ot  tltl 0.  body  o t  abstract  and  indexing  annotation  must  ftf  entered  whan  tha  overall  report  It  clas tilled) 

i  ORICINATINC  activity  (Corporate  author) 

I  The  Aerospace  Corporation 
(  El  Segundo,  California 

Ze  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

26  CROUP 

j  J  REPORT  TITLE 

IMAGING  OF  EXTENDED  OBJECTS  THROUGH  A  TURBULENT 

1  ATMOSPHERE 

1 

*  OE3CfliP"'iVE  note>  (Type  ot  tepott  end  Inclusive  date*) 

Richard  F.  j.  utomirski  and  Hal  T.  Yura 


6  REPORT  DATE 

7  3  TIIM  19 

m  ngggu 

76  NO  OF  REFS 

11 

9«  CONTRACT  OR  GRAN?  NO 

F04701-72-C-0073 

6  PROJECT  NO 

9a  ORIGINATOR'S  REPORT  numbc Rf S ) 

TR-0073(3230-30)-4 

c 

1  d 

96  OTHER  REPORT  NOfS;  (Any 
this  report) 

SAMSO-TR-73- 

other  numbers  that  may  be  assigned 

213 

•  0  DISTRIBUTION  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

II  SUPPLEMENTARY  NOTES 

«3  'ABSTRACT 

12  SPONSORING  MILITARY  ACTIVITY 

Space  and  Missile  Systems  Organization 

Air  Force  Systems  Command 

Los  Angeles,  California 

I  The  extended  -luygens -Fresnel  principle  is  used  to  derive  an  explicit 

expression  for  the  image  plane  illuminance  distribution  of  an  extended  object 
I  with  an  arbitrary  luminance  emittance  distribution.  The  combined  simul- 
I  taneous  effects  of  attenuation,  background  luminance,  and  atmospheric 
I  turbulence  are  given.  A  quantitative  comparison  of  these  effects  is  made, 
j  and  their  contribution  U  the  overall  loss  in  resolution  is  given.  In  particular, 
I  'We.  derive  a  quantitative  expression  for  the  point  spread  function  of  the  com- 
I  bined  atmospheric-optical  system.  Explicit  expressions  arc  derived  for  both 
?  the  atmospheric-optical  system  modulation  transfer  function  (MTF)  and 
f  image  plane  modulation.  Numerical  results  for  the  image  plane  modulation 
-  are  presented  for  imaging  both  up  and  down  along  an  atmospheric  path. 

I 

£ 

j 

I 

I 

il _ _ 

\<u 


(DD  F0RM  1473 

IFACSIMILCI 


UNCLASSIFIED 
Security  Classification 


UNCLASSIFIED 
Security  Classification 


Coherence 

Imaging 

Light  propagation 
Random  media 


UNCLASSIFIED 
Security  Classification 


Air  Force  Report  No. 
SAMSO-TR-73-213 


Aerospace  Report  No. 
TR-0073(3230-30)-4 


IMAGING  OF  EXTENDED  OBJECTS 
THROUGH  A  TURBULENT  ATMOSPHERE 


Prepared  by 
R.  F.  Lutomirski 

Pacific  Sierra  Research  Corporation 
and 

H.T.  Yura 

Electronics  Research  Laboratory 
The  Aerospace  Corporation 


73  JUN  28 


Laboratory  Operations 
THE  AEROSPACE  CORPORATION 


Prepared  for 

SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
AIR  FORCE  SYSTEMS  COMMAND 
LOS  ANGELES  AIR  FORCE  STATION 
Los  Angeles,  California 


Approved  for  public  release; 
distribution  unlimited 


FOREWORD 


This  report  is  published  by  The  Aerospace  Corporation,  El  Segundo, 
California,  under  Air  Force  Contract  No.  F04701-72-C-0073. 

This  report,  which  documents  research  carried  out  from  October 
through  December  1972,  was  submitted  16  April  1973  to  Lieutenant  Ernest 
L.  Lockwood,  DYAE,  for  review  and  approval. 


A.  H.  Silver,  Director 
Electronics  Research  Laboratory 


Publication  of  this  report  does  not  constitute  Air  Force  approval  of 
the  report's  findings  or  conclusions.  It  is  published  only  for  the  exchange 
and  stimulation  of  ideas. 


Ernest  L.  Lockwood 
IstLt.  ,  U.S.  Air  Force 
Project  Officer 


ABSTRACT 

The  extended  Huygens -Fresnel  principle  is'  used  to  derive  an  explicit 
expression  for  the  image  plane  illuminance  distribution  of  an  extended  object 
with  an  arbitrary  lum inane c  emittance  distribution.  The  cpmbined  simul¬ 
taneous  effects  of  attenuation,  background  luminance,  and  atmospheric  tur¬ 
bulence  are  given.  A  quantitative  comparison  of  these  effects  is  made,  and 
their  contribution  to  the  overall  loss  in  resolution  is  given.  In  particular, 
we  derive  a  quantitative  expression  for  the  point  spread  function  of  the 
combined  atmospheric-optical  system.  Explicit  expressions  are  derived 
for  both  the  atmospheric-optical  system  modulation  transfer  function  (MTi) 
and  image  plane  modulation.  Numerical  results  for  the  image  plane  modula¬ 
tion  are  presented  for  imaging  both  up  and  down  along  an  atmospheric  path 
under  various  viewing  conditions  . 
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I.  INTRODUCTION 


The  effects  of  the  atmosphere  significantly  degrade  the  imaging 
capability  of  optical  systems  in  many  applications.  The  mechanisms  that 
produce  image  degradation  include  the  attenuation  of  the  light  from  the  object 
of  interest,  light  scattered  into  the  imaging  system  from  aerosols  and 
unwanted  objects  (i.e. ,  background  luminance),  and  image  blurring  due  to 
atmospheric  turbulence.  The  limitations  on  imaging  capability  resulting 
from  each  of  these  mechanisms  acting  separately  has  been  treated 
previously.^  However,  to  the  best  of  our  knowledge  no  derivation  of  the 
image  plane  illuminance  of  an  extended  object,  which  includes  the  combined, 
simultaneous  degrading  effects  of  the  above  mechanisms,  exists  in  the 
literature. 

We  present  a  systematic  derivation  of  image  plane  illuminance 

and  the  corresponding  limits  of  resolution,  using  the  extended  Huygens - 

(2) 

F.resnel  principle.  This  technique  permits  a  quantitative  comparison  of 
the  contribution  of  each  of  the  above  mechanisms  to  the  overall  loss  of 
resolution.  In  particular,  an  explicit  expression  is  obtained  for  the  point 
spread  function,  defined  as  the  image  plane  illuminance  due  to  a  point  source 
located  in  the  object  plane,  of  the  combined  atmospheric-optical  system. 

The  extension  of  the  vacuum  Huygens -Fresnel  principle  to  the 
atmosphere  is  that  the  optical-field  at  an  arbitrary  point  in  space  can  again  be 
represented  as  a  superposition  of  spherical  waves  emitted  from  the  primary 
wavefront  existing  at  the  object,  but  each  spherical  wavelet  is  now  deter¬ 
mined  by  the  manner  in  which  it  propagates  through  the  medium.  As  a  result, 
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a  physical  interpretation  of  the  solution,  in  any  given  specific  case,  is 
easily  obtained.  The  use  of  this  principle  gives  directly  the  solutions  to 
many  problems  of  interest.  ^.s  an  additional  illustration  of  the  usc- 

fulness-of  this  principle,  we  apply  it  here  to  the  calculation  of  image  for¬ 
mation  through  the  atmosphere.  The  extended  Huygens -Fresnel  principle 
is  particularly  suited  for  this  application,  since  it  permits  the  separation 
of  the  geometry  of  the  problem,  i.e. ,  the  extended  source  distribution, 
from  the  propagation  problem,  which  is  determined  completely  by  the  man¬ 
ner  in  which  spherical  waves  propagate  from  various  elements  of  the  object 
through  the  atmosphere  to  the  imaging  optics. 

In  Sec.  II,  the  extended  Huygens-Fresnel  principle  is  used  to  derive 
an  explicit  expression  for  the  image  plane  illuminance  of  an  extended  object, 
which  includes  all  of  the  degradations  given  earlier.  In  Sec.  HI,  an 
expression  for  the  optical  transfer  functions  (OTF)  of  the  atmospheric- 
optical  system  is  obtained.  In  particular,  we  derive  explicit  expressions 
for  both  the  atmospheric-optical  system  modulation  transfer  function  (MTF) 
and  image  plane  modulation.  Finally,  in  Sec.  IV,  we  present  numerical 
results  for  the  image  plane  modulation  for  imaging  both  up  and  down  along 
an  atmospheric  path  under  various  viewing  conditions. 
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H.  GENERAL  CONSIDERATIONS:  IMAGE  PLANE  ILLUMINANCE 


We  consider  tlic  problem  of  imaging  an  extended  object  with  an 
arbitrary  luminance  cmittance  distribution  through  a  turbulent  atmosphere. 
Under  natural  illumination  conditions,  the  object  can  be  considered  to  be  an 
incoherent  radiator;  i-c.,  the  light  emitted  cither  by  reflection  or  direct 
emission  from  each  point  of  the  object  is  uncorrclatcd  with  that  emitted  from 
any  other  point.  Even  when  an  object  is  illuminated  by  laser  light,  the 
reflected  light  is  emitted  as  essentially  incoherent  radiation  when  the  scale 
of  irregularities  of  the  object  is  small  compared  with  the  coherence  length  of 
the  light  beam.  (This  is  the  case  for  almost  all  objects  of  practical  interest. ) 

In  the  configuration  depicted  in  Fig.  1,  the  object  is  located  a  dis¬ 
tance  z  from  a  thin  positive  lens  of  focal  length  f,  and  we  want  to  determine 
the  resulting  image  illuminance  distribution  in  the  image  plane  of  the  lens. 
The  distance  f'  of  the  image  plane  from  the  lens  plane  is  given,  for  a  focused 
system,  by 


The  turbulent  medium  between  the  object  and  lens  is  characterized 
by  the  spectral  density  of  the  index  of  refraction  fluctuations  ^  and  a 
wavelength-dependent  attenuation  constant  a,  which  is  the  sum  of  the  absorp¬ 
tion  coefficient  and  a  Rayleigh,  Mie,  and  nonselective  scattering  coeffi- 
(7) 

cient. 1  To  account  for  this  attenuation,  quantities  proportional  to  the 
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TURBULENT 

ATMOSPHERE 


square  of  the  field  contain  the  transmission  factor  exp(-crz).  In  addition  to 
the  light  emanating  from  the  object,  we  assume  that  the  lens  is  illuminated 
by  some  background  light  that  is  incoherent  with  respect  to  the  light  coming 
from  the  object.  The  resulting  optical  illuminance  entering  the  lens  is  then 
the  sum  of  that  coming  from  the  object  and  the  background  light. 

In  addition,  we  assume,  for  simplicity,  a  two-dimensional  object, 
which  is  described  by  a  given  field  distribution  UQ(v),  where  v  is  an  object 
coordinate.  Coordinates  in  the  lens  and  image  plane  will  be  denoted  by  the 
two-dimensional  vectors  r  and  x,  respectively. 

In  the  following,  the  resulting  image  illuminance  U(r)  at  the  entrance 
pupil  of  the  lens  is  obtained  from  a  direct  application  of  the  extended  Huygens- 
Fresnel  principle.  The  modification  of  the  results  of  Ref.  2  to  include  atmo¬ 
spheric  attenuation  is^ 

U(r)  =  exp  (-  ^)(- %)fciZ,r)  Uo(v)  d2y  (2) 

where  the  integration  extends  over  the  object  distribution,  G  is  the  Green's 
function 


G(v,  r) s  (~)  exp[iks(v,  r)  +  'My.,!)]  (3) 

s(y,_r)  is  the  geometrical  distance  between  the  object  pointy  and  the  lens 

(2) 

point  _r,  >|<  is  the  complex  phase  perturbation  due  to  the  turbulent  medium, 
and  k  is  the  optical  wave  number.  We  take  the  spatial  extent  of  the  object 
(and  lens)  to  be  much  less  than  the  propagation  distance  z,  from  which  it 
follows  that 


(4) 


s(v,r) 


z  f 


(v  -  r)2 

2  Z 


The  field  V  (r)  at  the  exit  surface  of  the  lens  is  given  by  V  (r)  = 

W(r)  exp[-ikr^/2f]  U(r).  The  quantity  W(r),  the  lens  factor,  is  unity  for 
jr|  <  D/2,  where  D  is  the  diameter  of  the  lens,  and  zero  otherwise,  and  the 
exponential  factor  is  just  the  transmission  factor  of  a  positive  lens  of  focal 
length  f.  The  instantaneous  field  distribution  in  the  image  plane  is  given  by 
the  vacuum  Huygens-Fresnel  principle 

Uj(x)  =  (-^)/v(r)Go(r,x)d2r 

=  (-§)/w(l)  exp[-i||-]u(r)d2r  (5a) 

and  substituting  Eq.  (2)  into  Eq.  (5a)  yields 

Uj(x)  =  (-||)  oxp(^-~jJd\ f  j2£  u0W  g<£-£)  G0(r,x) 

X  W  (r )  exp^-^j  (5b) 

where  the  r  integration  extends  over  the  lens  and 

G0(£'£>  s(f)  expfiksfr,  x)] 

The  instantaneous  illuminance  distribution  in  the  image  plane  is 

given  by 
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Ei(x)  =  |U.(x}|2 

=  (-%)  exp(-orz) J <lZv  J d2.r J d2r'  UQ(i)  U*(£!)  G(jr,r) 

2  2 

X  cljr’.r1)  Go(r,x)  G*(r',x)  W(r)  W(r')  exp [- 

Since  the  object  radiates  incoherently,  the  only  contribution  to  the 
integrals  in  the  above  are  obtained  for  £  ~  yr*m  AVe  have 

<E.(x)>  =  Aexp(-orz)  J d2jr  J ^ L  J "d*V  Eq(x)  <G(x,r)  G'(jr,  r_' )> GQ(_r,x)  G^r'.x) 
X  W(r)  W(r')  exp[-  (?> 

2  2 

where  EQ(y)  =  |UQ(y)|  is  the  luminous  omittance  (lm/m  )  distribution  of  the 
object,  and  since  we  are  dealing  with  a  random  medium,  an  ensemble  aver¬ 
age  (denoted  by  angular  brackets)  has  been  taken.  In  the  appendix,  we  show 

that  by  considering  the  vacuum  limit  of  the  illuminance  when  the  object  is 

2  3 

reduced  to  a  point,  the  constant  can  be  determined  as  A  =  k  /8n  . 

Setting  G  =  Gq  expW>),  we  obtain 

<G(i,  r.)  G*(x>  r')>  =  GQ(x,r)  G^(x,  r')  <exp[if>(x.,  r)  +  «|»*(£,r')]>  (8) 

where 

G0(x,  r)  =(^-)exp[iks(x,r)]  (9) 
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The  quantity  <exp[ii(£,  rj  +  4,*(^.£.,>]^  is  recognized  as  the  mutual ' 

coherence  function  (MCF)  of  a  spherical  wave  propagating  in  the  turbulent 

(2) 

medium  from  the  point  on  the  object  to  the  lens  points  r_  and  r).  We 
assume  here  that  the  MCF  is  independent  of  source  coordinate  y_;  that  is, 
we  assume  the  region  of  isoplanatism  is  larger  than  the  spatial  extent  of 
the  object.  Further,  for  the  case  of  homogeneous  isotropic  turbulence  the 
MCF  can  be  expressed  as  Mg(p),  where  p  =|r^  -  _r"| .  The  turbulence  prop¬ 
erties  of  the  medium  are  contained  in  the  spherical  wave  MCF,  and  thus  a 
knowledge  of  Mg  is  sufficient  for  the  determination  of  the  resulting  average 
image  illuminance  characteristics. 

Substituting  Eqs.  (6),  (8),  and  (9)  into  Eq.  (7)  and  changing  vari¬ 
ables  from  jr  and  r)  to  £  =  r_  -  r_'  and  R.  =  (£  +  r))/Z,  the  average  illuminance 
can  be  written  in  the  form 

<E.(x)>  =  J S(£  +  xz/f)  Eq(z)  (10a) 


where  the  point  source  response  function  S  is  given  by 

Six)  =  — exP(-ttZ> /  ml(d)  ms(p>  exP  [-  lkz£,J]d2g.  <10b> 

and  M^  is  the  lens  transfer  function  [Eq.  (17)]. 

Equations  (10a)  and  (10b)  reveal  that  the  average  illuminance  can  be 
expressed  as  a  convolution  of  the  point-source  response  function  with  the 
luminous  emittance  distribution  of  the  object,  the  point  source  response 
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function  is  proportional  to  the  Fourier-transform  of  the  product  of  the  lens 
transfer  function  and  the  spherical  wave  MCF. 

In  addition  to  the  illuminance  distribution  resulting  from  the  object, 
we  must  add  the  corresponding  illuminance  distribution  due  to  the  background 

light.  Background  light  can-be  characterized  by  a  uniform  luminance 
2 

B^(lm/m  -  sr),  resulting  in  an  image  plane  illuminance  E^. 


Eb  =  Bb 


(area  of  lens) 


(distance  lens  to  image  plane) 


Bb(irD  /4)  irBb 


(f1) 


4F 


(li) 


where  the  F  number  of  the  lens  is  the  ratio  of  the  focal  length  to  the  lens 
diameter.  The  background  image  brightness  is  uniform  over  the  image 
plane  and  is  independent  of  the  distance  to  the  object.  The  F  number  of  the 
lens  determines  the  image  brightness;  the  smaller  the  F  number,  the  brighter 
the  image.  From  Eq.  (10b),  this  applies  to  object  light  as  well  as  the  back¬ 
ground  light. 

Thus,  the  total  illuminance  distribution  in  the  image  plane  is 

Et(x)  =  <Ej(x)>  +  Eb  (12) 


Background  luminance  (or  path  luminance)  depends  on  the  relative 
positions  of  the  sun,  the  observer,  and  the  scene.  It  also  depends  on  the 
observer's  altitude  and  the  condition  of  the  atmosphere;  it  is  different  for 


an  observer  on  the  ground  looking  at  an  airplane  than  for  the  reverse 
situation.  Figure  2  is  a  schematic  of  the  viewing  geometry.  The  param¬ 
eters  indicated  in  the  figure  are:  observer's  altitude  H,  zenith  angle  between 
the  local  vertical  and  the  light  of  sight  Q,  azimuth  angle  between  the  line  of 
sight  and  the  observer  sun  line  measured  in  a  horizontal  plane  <p,  and  solar 
zenith  angle  p.  Standard  or  average  values  of  the  path  luminance  have  been 
tabulated  for  a  solar  zenith  angle  of  P  =  41. 5  deg  and, an  "average"  atmos- 

/QV 

phere.  '  Figure  3  is  adapted  from  tables  in  Ref.  8. 
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Fig.  3a.  Average  path  luminance  as  a  function  of  altitude  and 
viewing  angle  for  p  =  41.5  deg  and  <p  =  0  deg 
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ZENITH  ANGLE  0  BETWEEN  OBSERVER  AND  TARGET  (deg) 


Fig.  3b.  Average  path  luminance  as  a  function  of  altitude  and 
viewing  angle  for  f}  =  41. 5  deg  and  6  =  90  deg 
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Fig.  3c.  Average  path  luminance  as  a  function  of  altitude  and 
viewing  angle  for  0  =  41. 5  deg  and  $  =  180  deg 
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in.  OPTICAL  TRANSFER  FUNCTION 

A  quantity  of  particular  interest  is  the  optical  transfer  function 
(OTF)  of  the  system.  This  quantity  is  defined  as  the  complex  (image)  illumi¬ 
nance  distribution  corresponding  to  a  unit  amplitude  (object)  luminous  omit¬ 
tance  of  a  single  spatial  frequency.  The  OTF  describes  how  both  the  ampli¬ 
tude  and  phase  of  the  image  distribution  change  as  a  function  of  spatial  fre¬ 
quency.  The  resolution  of  the  viewing  system  is  determined  directly  from 
the  OTF.  When  only  amplitude  changes  arc  considered,  as  is  often  the  case, 
the  OTF  is  equal  to  the  MTF. 

In  this  regard,  v.-c  arc  led  to  consider  the  two-dimensional  Fourier 
transform  of  the  image  illuminance  distribution.  The  illuminance  consists 
of  two  terms,  the  background  term  and  that  due  to  light  from  the  object. 

Since  the  background  term  E,  is  constant,  its  Fourier  transform  is  just  a 
dc  term  at  zerq  spatial  frequency.  The  Fourier  transform  of  is 


where 


£. 


(15) 


and 

Eq(Q)  = f  Eo(v)  exp(-iQ-i)  d\  (16) 

is  the  Fourier  transform  of  the  object  distribution.  The  minus  sign  appear¬ 
ing  in  Eq  in  Eq.  (14)  corresponds  to  the  inversion  of  the  image  by  the  lens 
while  the  factor  f'/z  corresponds  to  the  linear  magnification  of  the  spatial 
distribution  of  the  object.  For  cases  of  interest  here,  z  is  large  compared 
to  the  focal  length  f,  and  “•■is  V  as  f. 

The  integral  in  the  bracket  in  Eq.  (14)  is  the  convolution  of  the  lens 
pupil  function  with  itself  and  is  given  by 

(^)  /dZ«.  W(H  f  I)  W(H  - 1)  =  ml(d)  <17> 

where  M^»  the  MTF  of  the  circular  lens,  ^is^  given  by 

Ml(d)  =  I  jCOS'1(D)-(D)[1  -(d)2]  }  (18) 


and  p  is  given  by  Eq.  (15).  In  the  general  case,  as  defined  by  Eq.  (17), 

is  the  MTF  of  a  lens  specified  by  an  arbitrary  pupil  function  W(rj. 
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The  OTF  is  given  by  the  ratio  of  the  object  to  image  Fourier 
components.  In  this  case,  the  OTF  is  real,  and  is  identical  to  the  MTF. 
From  Eqs.  (14)  and  (17),  we  find 

MTF  =  faW  exp(‘"z)  Ms(p)  Ml(d)  (19) 

where  p  cr  Kf/k  =  K'f\  where  K1  =  K/2tt  has  the  dimensions  cycles  per  unit 
length.  In  this  case,  the  MTF  of  the  system  is  the  product  of  the  spherical 
wave  MCF  of  the  medium,  evaluated  at  p  =  Kf/k,  and  the  MTF  of  the  lens, 
both  of  which  are  functions  of  spatial  frequency.  In  addition,  the  MTF 
contains  the  factor  exp(-az),  which  is  independent  of  spatial  frequency  and 
corresponds  to  the  loss  of  signal  illuminance  due  to  large  angle  scattering 
and  absorption.  The  decomposition  of  the  overall  system  MTF  into  a 
product  of  independent  terms  (e.g. ,  medium,  lens)  is  possible  only  if  the 
spherical  wave  MTF  of  the  medium  is  independent  of  the  origin  of  the 
spherical  wave  source,  i.  e. ,  the  entire  object  lies  in  an  isoplanatic  patch. 

As  an  illustrative  example,  we  calculate  the  resulting  degradation 
of  contrast  of  a  sinusoidal  object.  We  consider  an  object  that  consists  of 
a  series  of  black  and  white  linps  or  some  other  repetitive  luminance 
pattern  of  a  given  spatial  fijepuency.  For  mathematical  simplicity,  we 
model  such  a  pattern  by  the  fpllowing  one-dimensional  function 


E 


E0(y)  =  (l  +  CRS  Ky] 


(20) 


(24) 


where 


4ttB. 
_ b 

n  '  Eq  oxp(-ffz) 


(25) 


Now,  for  a  Lambertian  surface  of  average  reflectivity  r,  the  average 
luminance  Eq  of  the  object  is  given  by 


(26) 


where  Eg  is  the  incident  illuminance.  By  substituting  Eq.  (26)  into 
Eq.  (25),  we  obtain 


4  ~  Egr  exp(-orz) 


(27) 


The  dimension  loss  quantity  q  is  independent  of  spatial  frequency  and  F 
number.  The  factor  (1  +  q)"*  represents  an  overall  reduction  in  image  mod¬ 
ulation.  Only  for  the  case  q  «  1  is  the  image  modulation  given  by  the  product 
of  the  spherical  wave  and  lens  MTF. 


IV.  NUMERICAL  EXAMPLES:  SYSTEM  RESOLUTION 

As  an  illustrative  example,  we  consider  the  case  of  looking  both  up 
or  down  with  a  6-in.  lens  from  an  altitude  of  3  km  along  a  slant  path  of  10  km. 
To  be  definite,  we  consider  viewing  a  solar  illuminated  horizontal  object  in 
the  visible.  We  set  <p,  the  azimuth  angle  between  the  line  of  sight  and  the 
observer  sun  line  measured  in  a  horizontal  plane,  equal  to  90  deg,  and  p, 

the  solar  zenith  angle,  equal  to  41.5  deg.  From  Fig.  4,  we  find  that  for 

3  2 

looking  down  (6  =  108  deg)  B^  =  3  X  10  cd/m  and  for  looking  up  (0  =  72 
3  2 

deg)  =  2  X  10  cd/m  .  For  this  viewing  geometry,  we  find  from  Ref.  9 
that  Eg  s  7.44  X  10^  lm/m^,  and  assuming  an  average  reflectivity  of  0.  5 
we  obtain,  from  Eq.  (27)  and  Ref.  9,  Figs.  4a  and  4b  and  the  Table  I. 

In  Figs.  (4a)  and  (4b),  we  have  plotted  the  image  modulation  look¬ 
ing  down  and  looking  up  as  a  function  of  angular  spatial  frequency  K'f  (which 
is  equal  to  p/\).  In  obtaining  these  curves,  we  have  used  Coulman's  data^^ 
for  the  temperature  structure  constant  profile  for  average  conditions  during 
■the  midpart  of  sunny  days , 


Table  I.  *1  For  Various  Atmospheric  Conditions 


Looking 

Extra  Clear 

Very  Clear 

■mi 

Down 

0,  92 

1.40 

6.60 

Up 

0.62 

0.93 

4.40 

Preceding  page  blank 
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For  any  given  receiving  system,  there  is  some  threshold  modulation 
mQ  (<1),  such  that  if  the  image  modulation  is  less  than  m0,  resolution  is  not 
possible.  Generally,  mQ  is  a  function  of  spatial  frequency  and  level  of 
illumination  reaching  the  receiver.  Here,  for  illustrative  purposes,  we 
arbitrarily  assume  mQ  is  a  constant  and  equal  to  5%.  In  this  case,  we  see, 
from  Fig.  4a,  that  for  a  very  clear  atmosphere  the  limiting  angular  resolu¬ 
tion  is  ~0.5  prad.  Conversely,  for  looking  up  under  these  same  conditions, 
the  limiting  resolution  is  ~3  prad.  (Both  of  these  limits  apply  to  resolution 
in  a  direction  normal  to  the  line  of  sight. ) 
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Fig.  4a.  Image  plane  modulation  for  summer  daytime 
turbulence- conditions,  looking  up 
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APPENDIX:  EVALUATION  OF  THE  CONSTANT  A 


To  evaluate  the  constant  A  in  Eq.  (7),  we  consider  the  vacuum 
limit  (i.  e. ,  or  =  0,  G  =  1)  of  the  image  plane  illuminance  as  the  diameter 
of  the  object  Ay  tends  toward  zero  (positioned  at  y  =  0)  and  of  luminous 
strength  tends  to  W„;  i.e. , 


W0  S  Ay— Q /Eo<i)  d\ 


In  this  case,  Eq.  (7)  becomes 


E.(x)  =  ~  | / Go(r,x)  W(r)  exp^^^y'--^— ']  d2r 


The  vacuum  propagator  Go(:r,x)  is  given,  in  the  paraxial  approximation,  by 


G0<E.x)  S  (jr)  expjikf  +  If'  ^  ]| 


Hence,  from  Eqs,  (A-2)  and  (A-3),\ve  obtain 


E.(x)  =■ 


-l|/exP(-^)w(r)d2r 

SI2"  /J  o(~r)MZ 


(Zf<ri  > 

?  2 

AW0(ira  )  f2J1(kax/fl) 
/,m2  (kaoc/f) 


Preceding  page  blank 


j 


where  a  is  the  radius  of  the  lens. 

The  total  luminous  flux  W,  impinging  on  the  image  plane  is  obtained 
by  integrating  Eq.  (A-4)  over  all  x.  We  obtain 


w,  =y E.(x)  d2x 


=  Zit  I  E.(x)  xdx 
*/0  ' 


ir3  AWq  D2 

kV" 


(A-5) 


where  D  =  2a  is  the  diameter  of  the  lens. 

The  quantity  Wj  is  equal  to  the  total  flux  impinging  on  the  lens, 
which  for  a  point  source  radiating  into  a  half  space  must  also  equal  the  source 
strength  WQ  multiplied  by  the  fractional  area  intercepted  by  the  lens  on  a 
hemisphere  of  radius  z: 


By  comparing  Eqs.  (A-5)  and  (A-6),  we  obtain 


(A-6) 


(A-7) 
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